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What gravitational field is generated by a massive quantum system in a spatial superposition? 

This is one of the most important questions in modern physics, and after decades of intensive the¬ 
oretical and experimental research, we still do not know the answer. On the experimental side, the 
difficulty lies in the fact that gravity is weak and requires large masses to be detectable. But for 
large masses, it becomes increasingly difficult to generate spatial quantum superpositions, which 
live sufficiently long to be detected. A delicate balance between opposite quantum and gravitational 
demands is needed. Here we show that this can be achieved in an optomechanics scenario. We 
propose an experimental setup, which allows to decide whether the gravitational field generated by 
a quantum system in a spatial superposition is the superposition of the two alternatives, or not. 

We estimate the magnitude of the effect and show that it offers good perspectives for observabil¬ 
ity. Performing the experiment will mark a breakthrough in our understanding of the relationship 
between gravity and quantum theory. 


Quantum field theory is one of the most successful the¬ 
ories ever formulated. All matter fields, together with the 
electromagnetic and nuclear forces, have been success¬ 
fully embodied in the quantum framework. They form 
the much celebrated standard model of elementary par¬ 
ticles, which not only has been confirmed in all advanced 
accelerator facilities, but has also become an essential 
ingredient for the description of the universe and its evo¬ 
lution. 

In light of this, it becomes obvious to seek a quantum 
formulation of gravity as well. Yet, the straightforward 
procedure for promoting the classical field as described by 
general relativity, into a quantum field, does not work. 
Over the decades, several strategies have been put for¬ 
ward, which turned into very sophisticated theories of 
gravity, perhaps the most advanced being string theory 
and loop quantum gravity. Yet, none of them has reached 
the goal of providing a fully consistent quantum theory 
of gravity. 

At this point, one might wonder whether the very idea 
of quantizing gravity is ill-posed mm- At the end of the 
day, according to general relativity, gravity is rather dif¬ 
ferent from all other forces. Actually, it is not a force at 
all, but a mere manifestation of the curvature of space- 
time, and there is no obvious reason why the standard 
approach to the quantisation of fields should work for 
spacetime as well. A future unified theory of quantum 
and gravitational phenomena might require a radical re¬ 
vision not only of our notions of space and time, but 
also of (quantum) matter. This scenario is growing in 
likeliness. 

From the experimental point of view, it has now been 
ascertained that quantum matter (i.e. matter in a gen¬ 
uine quantum state, such as a coherent superposition 
state) couples to the Earth’s gravity in the most obvious 


way. This has been confirmed in neutron [3], atom [4] 
interferometers and used for velocity selection in molec¬ 
ular interferometry [^. However, in all cases, the grav¬ 
itational field is classical, i.e. it is generated by a dis¬ 
tribution of matter (the Earth) in a fully classical state. 
Therefore, the plethora of successful experiments men¬ 
tioned above does not provide hints, unfortunately, on 
whether gravity is quantum or not. 

The large attention and media coverage about the 
BICEP2 experiment having shown the quantum origin 
of primordial gravitational fluctuations [B], subsequently 
disproved by Planck’s data analysis [7], testifies the im¬ 
portance and urgency of a pragmatic assessment of the 
question of whether gravity is quantum or not. 

In this paper, we propose an approach where a meso¬ 
scopic system is forced in the superposition of two dif¬ 
ferent positions in space, and its gravitational field is 
explored by a probe (Fig. 1). Using the exquisite poten¬ 
tial for transduction officered by optomechanics, we can 
in principle determine whether the gravitational field is 
the superposition of the two gravitational fields associ¬ 
ated to the two different states of the system, or not. 
The first case amounts to a quantum behavior of grav¬ 
ity, the second to a classical-like one. We show that the 
sensitivity necessary to appreciate the difference between 
such behaviors is close to the current state of the art in 
specific optomechanical configuration, although quite de¬ 
manding. 

Framework.- We consider a system SI (with mass mi) 
prepared in a superposition of two different positions in 
space. The wave function is V’(ri) = («(ri) + ,d(ri)) 

with 7(ri) = (rijy) (7 = a, 13) and {a\(3) = 0, stat¬ 
ing the distinguishability (in a macroscopic sense) of the 
two states. SI generates a gravitational field that can be 
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FIG. 1: A two-body system that interacts only gravitation¬ 
ally. To study how these systems interact gravitationally, 
we only need to consider the degrees of freedom in a two- 
dimensional xy-plane. System SI is initially prepared in a 
superposition of two wavepackets delocalised along the x axis 
of the reference frame (red balls). System S2 (the probing sys¬ 
tem) is initially prepared in a narrowly localized wavepacket 
(blue ball). The parameters of the experiment are adjusted 
in such a way that the non-gravitational interactions between 
SI and S2 are, for all practical purposes, negligible. 


probed with the help of a second system S2 (with mass 
m 2 )- The latter is prepared in a localized state with asso¬ 
ciated wave function (/)(r 2 ) [cf. Fig. 01- The question we 
address is: which gravitational field does S2 experience, 
as generated by SI being in a spatial superposition? We 
consider two alternatives. 

Quantum gravity scenario.- Although we do not yet have 
a fully consistent and complete theory of quantum grav¬ 
ity, we can safely claim that, should gravity be quantum, 
this would be manifested in SI generating the superpo¬ 
sition of two weak gravitational fields, one associated to 
the state Q!(ri) and the other one to the state /3(ri). Lin¬ 
earity is the characteristic trait of quantum theory, which 
should find place in any quantization recipe for gravity. 
On top of this, nonlinear quantum theories very easily 
run into serious troubles [M]. Therefore SI generating 
the superposition of two gravitational fields is the most 
likely scenario in a world where also gravity is quantum. 

Then S2 reacts to the situation described here above by 
turning into a superposition of being attracted towards 
the region A where Q;(ri) is different from zero, and be¬ 
ing attracted towards B, where /3(ri) is different from 
zero. If the motion of S2 is constrained within an har¬ 
monic trap, a position measurement will reveal it being 
displaced slightly away from equilibrium, alternatively 
closer to A or B. 

Semiclassical gravity scenario.- Suppose instead that 
gravity, for whatever reason, is fundamentally classical. 
In this case, no one really knows which gravitational 
field SI generates. However, a natural answer is that 
the (square modulus of the) wave function acts as some 


sort of matter distribution, which generates the gravita¬ 
tional field. This is what semi-classical Einstein equa¬ 
tions [in [H [M] predict, if taken seriously as funda¬ 
mental equations. It is also compatible with our current 
knowledge of gravity. 

In this case, the gravitational field is not the quan¬ 
tum superposition of the two quantum fields associated 
to Q:(ri) and /3(ri), but the classical sum of the two classi¬ 
cal fields generated by their square modulus, respectively. 
In such conditions, S2 feels a force which pulls it in be¬ 
tween region A and B. Quite evidently, the two cases 
imply two different motions for the probe S2, such differ¬ 
ence being the way to discriminate between classical-like 
and quantum treatment. 

Theoretical modelling.- We refer again to the situation il¬ 
lustrated in FigEl In what follows, when no explicit time 
dependence is reported, we imply t — 0. We let the two 
systems interact gravitationally for a time r, and then 
measure the position of S2 along the x-axis. The time- 
scale T can be, at most, the lifetime of the superposition 
state tp. Also, experimental parameters are adjusted such 
that all interactions, except gravity, are negligible, for all 
practical purposes. 

In the quantum approach, the total Hamiltonian is 
given by iL = Hi + H 2 + HL 2 (ri - r 2 ), where Hi ,2 
are the Hamiltonians of SI and S2, respectively, and 
Fi 2 (ri —r 2 ) is the Newtonian interaction. The final 
state of the overall system is given by 'I'(ri,r 2 ,r) = 
[tI'a(ri,r 2 ,T)-|-'I',s(ri,r 2 ,T)]/ , where each term is the 

solution of the equation ihdt'^a,p = with ini¬ 
tial conditions r 2 ) = a(ri)(/)(r 2 ) and 'I'^(ri,r 2 ) = 

/3(ri)(^(r2). We assume mi ^ m 2 , implying an adiabatic 
approximation in which the degrees of freedom of the two 
systems can be separated as 


T'.^(ri,r2,t) = 7(ri,t)(^^(r2,t), (7 = 0,^), (1) 


where the motion of a(ri,t) and /I(ri,t) are determined 
by Hi, while (j)j{r 2 ,t) evolves with the Hamiltonian H^ = 
H 2 + Vy with 

Vy = -Gmim2 (7 = 0 ,/?). ( 2 ) 

J |ri-r2| 


In this quantum scenario, the initial superposition state 
of SI generates a superposition of gravitational helds, 
which in turn generates a superposition of motions for 
S2. 

On the other hand, in the semiclassical treatment of 
gravity, and under the same approximations discussed 
above, the evolution of S2 is determined by H^, = H 2 + 
V(,i, where the gravitational potential now reads 


Ki = —Gmi m 2 




ri 


|ri -r 2 | ■ 


( 3 ) 


Here, the evolution of ip{ri,t) is determined by Hi. 
Eq. (§ shows the signatures of a classical treatment of 
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FIG. 2: (a) The gravitational field acting on S2 is a linear 
combination of gravitational fields produced by SI being in a 
superposeed state, (b) The semi-classical treatment of grav¬ 
ity, where the gravitational field acting on S2 is the one pro¬ 
duced by a total mass mi with density | (|Q:(a;)|^ + |/3(a;)|^). 




gravity: the wave function gives the mass density, and 
there is no quantum superposition of gravitational fields. 

One can further approximate the above mentioned 
gravitational potentials, which will be quite useful when 
we compute the motion of optomechanical systems. 
Henceforth we will work in the Heisenberg picture. We 
assume that the quantum fluctuations around the mean 
values for SI are small. Therefore, Vy in Eq. (|^ can be 
approximated as 


A ^ Gmim 2 


(7 = a,/ 3 ). 


( 4 ) 


Assuming that the quantum fluctuations around the 
mean values for S 2 are also small, V-y can be expanded in 
Taylor series as 




Gmim 2 


|(ri(t))^ - (r2(t))0| 


e Gmim2((ri(f))^ - {r2it))^) 


( 5 ) 


where Sr2{t) = r2{t) - {r2{t))^. 

The same procedure can applied to Ki- As for the 
quantum case, assuming that the fluctuations in the mo¬ 
tions of SI and S 2 are small, we find 


K, 


- E 

+Sr 2 (t) ■ 


Gmim 2 

2|(ri(t))^ - (i'2(t)),^| 

GTOiTO2((ri(t))^ - {r2{t))<p) 
2|(ri(t)),-(r2(t))^|3 


( 6 ) 


Optomechanical test- We now consider the exquisite po¬ 
tential for motional transduction offered by optomechan- 
ics and let system S 2 be the movable end-mirror of an 
optomechanical cavity. On the other hand, we shall 
not specify explicitly what is the chosen embodiment for 
SI, which could well be a second vibrating mechanical 
structure. Explicit configurations will be described else¬ 
where im. 

The transduction cavity is pumped by an external laser 
field and S 2 is in contact with a bath of phononic modes. 


FIG. 3: The proposed set-up for the optomechanical falsifi¬ 
cation of quantum/classical gravity. A system SI is prepared 
in a superposition of two localised states at ±dx along the 
X axis. An optomechanical cavity acts as transducer and a 
probe of (potentially quantum) gravity effects S2: the effect 
of the gravitational coupling between SI and the mechanical 
oscillator of an optomechanical cavity induces an effect on 
the variance of the position fluctuations of the oscillator. The 
mean position of the latter along the x axis is X 2 - The cavity 
is pumped by an external field (frequency cjq and coupling 
rate £). 


The axis of the optomechanical cavity, as well as the de- 
localization axis of SI, is assumed to lay along the x-axis 
of a reference frame. Denoting the position operator of 
S 2 along the cavity axis by X2, and its momentum by 
P2, moving to a frame rotating at the frequency of the 
pumping field, we find the Hamiltonian model m 

H 2 = h{wc — oJo)a^CL — fixo)ax2 

1 2 2 P2 -kcf t ^ 

-I- - m2 Wo Xo + v-h ih£ (o' — a ), 

2 2m2 

where ujq is the frequency of the external laser, Wc is the 
frequency of the cavity mode derived by the laser, W2 
is the harmonic frequency of the mechanical oscillator, 
X = Wc/T is the optomechanical coupling constant be¬ 
tween the cavity and the mechanical oscillator with L 
the size of the cavity, and £ = 2 kV/ hojQ with V the 
laser power and k the cavity photon decay rate. Follow¬ 
ing conventional approach, we expand each operator as 
0 = 0 + SO with O the steady-state mean value and SO 
small quantum fluctuation around O. Accordingly, one 
finds: p2 = 0 , X2 = fi.x|dp/m2a;|, and a = £/{k + iA) 
with A = LUc — oJq — x^2- We now assume that the 
mean-value of the position of system S 2 , {r2{t))^, takes 
a steady-state value. This implies that the coordinates 
of S 2 in the Cartesian reference frame that we have 
chosen are {r2{t))^ « (x2,dy,0). Also, as mi S> m2, 
within the aforementioned adiabatic approximation we 
have (ri(t))a « (da;,0,0) and (r2(t))/3 » (-da:,0, 0) [cf. 
Fig. [^. We introduce these approximations into equa- 














































4 


tions of and V^i and, by taking dy ^ dx,X 2 , we find 




Gmim 2 


1 


X 2 


\ “y 

with j = a, /3, Sa = —sp = — 1 , and 


t Sjdx ^ If. 

^ ix. - - 5y, 


y,y 


Kl«- 


Gmim 2 


^ ^ Sy 2 

S dy 


( 8 ) 


( 9 ) 


The derivative of the potentials above with respect to 6 x 2 
contributes a term in the quantum Langevin equation 
of the momentum. Accordingly, the quantum Langevin 
equations read 

^-Sx 2 {t) = 6 p 2 {t)tm 2 , 


dt 

d 

dt 

d 

dt 


Sa{t) = —{iA + K) 6 a{t) + ixa Sx 2 {t) + Sai^{t), 


dp2{t) = —m2UJ2Sx2{t) + hxlaSa^ (t) +a*Sa{t)] 

- - lmSp 2 {t) + ^{t), 

( 10 ) 

where = 7 , cl. Notice that dV^fdSx 2 is not 
an operator-valued function. We shall denote / = 
—dVjd 5 x 2 - Solving the above equations in the frequency 
domain gives us 

5x2{uj) = ^{[A^ + {k- ia;)^][^(a;) -k 2 -KfJ{uj)\ 

+ihxV^ [d{iK + w—A)(5a)(,(a;) -I- d*{iK + a;-|-A)5ai„(a;)]} 

( 11 ) 

where D{uj) = m[A^ + (k — iujY'Wu)^ — wf -I- + 

2?ix^|ap A. The correlation functions of the noise oper¬ 
ators are 

(Sai„(uj) Sal^(ft)} = 2 tt6(uj + fl) 

(^(w) ^(n)) = 2'Kh'ym'muj [1 -I- coth(/ra;)] 5{(jJ -|- fl) 

with pi = h/ksT. All other correlators are zero. There¬ 
fore, one finds the spectrum of fluctuations in the position 
of the mechanical oscillator S2 as 

Sx 2 (w) = + K^+ 

-I- hmojXm coth {pLUj) [(A^ + + 4K^a;^]} 


( 12 ) 


-k 


D{oj)D{0) 


J(w)(A^ -|- [A^ -I- (k — iw)^]. 


( 13 ) 

The first two lines in this expression reproduce the stan¬ 
dard density noise spectrum of an optomechanical sys¬ 
tem. On the other hand, the term proportional to 
in Eq. (13) is the result of the gravitational interaction. 
This contribution can be directly observed in the variance 
of the position of S2. The variance of the fluctuations in 
the position of mechanical oscillator S2 is given by 

1 /'+°° 

(((5x2)^) = —^ dw5^3(u;). (14) 


Introducing Eq. (13) into the above expression yields 


{{ 6 x 2 )^) = {( 6 X 2 )% + 


D2(0) ’ 


(15) 


where ((i 5 x 2 )^)o denotes the variance of the position fluc¬ 
tuations of the mechanical oscillator when there is no 
gravitational interaction. We also used = —dVi,jddx 2 
where V can be either Ey, or E, [cf. Eqs. (|^ and ([^]. 
Explicitly 

dVj Gmim 2 . 9Ei Gmim 2 X 2 

^2"^ ^(^2 + s 7 dx); ^3 , 

(16) 


where 7 = q ;,/ 3 , and Sa = —sp = —1. As one can appre¬ 
ciate from Eq. (151, the gravitational interaction between 
SI and S2 manifests as an extra widening in the position 
distribution of S2. Eqs. (151 and (16) allow to evaluate 
the difference between a classical and a quantum treat¬ 
ment of gravity. As a figure of merit we can indeed take 


0 = Gmim 2 


(A^ -I- K^)^dx{dx + X2) 

dl\Dm 


(17) 


which is the difference between the standard deviation in 
classical and quantum cases, and has the dimension of 
a length. Our goal now is to achieve the largest possi¬ 
ble deviation. Upon inspection, one can see that 71(0) 
is minimized for A = 0 and uj = ui 2 - Moreover, by as¬ 
suming the (experimentally undemanding) sideband-not- 
resolved limit given by the condition k ^ a; 2 , we find the 
optimal expression 


0 * = Gmi 


\/ dx{dx + X 2 ) 
dluj 2 Xm 


Gmidx 
dluJ 2 Xm ’ 


(18) 


where the last expression is valid by assuming dx X 2 - 
This result shows that a high mechanical quality factor of 
a low-frequency oscillator would bring 0* to values close 
to observability, provided that the distance dy between 
the centres of mass of SI and S2 is larger than the linear 
dimension of the objects along the y axis. An estimate 
is as follows: we take £ = 6 x lO^^Hz, Xml"2-'K = lOOHz, 
AC = 9 X lO^Hz, a;c/27r = 3.7 x W2/27r = fO^Hz, 

mi = lOOng, and a cavity of 1mm length. Consis¬ 
tently with the formal approach above, we work under 
the assumptions dy dx,X 2 - A suitable range of val¬ 
ues for dy is from 10 “®m, which would be suited for 
micromechanical systems, to 10 “®m, which would imply 
the use of a nanomechanical oscillator (possibly embod¬ 
ied by a carbon nanotube or a graphene sheet, such as 
in Ref. [H]). Eor m 2 = Ing (m 2 = 3 x 10“^^g), we have 
X 2 « 7.7 X 10“^°m {x 2 ~ 2.6 x 10“^m). Eor dy = 10“®m 
{dy = 10“®m ) and = 5 x 10“^m {dx = 5 x 10“®m), 
we find 0* « 1.3 X 10“®G (0* « 1.3 x IQ-^G). The 
optimality of such value is assessed in Fig. where we 
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0/G [Kgs^/m^] 



FIG. 4: The red lines show the plot of 0 against the frequency 
uj for dy = 10“®m (lower red curve), and dy = 10“®m with 
m 2 = Ing (upper red curve). All other parameters as in the 
body of the paper. The horizontal blue lines show the values 
of 0*/G corresponding to our two choices of dy. We have 
maxi,; 0 = 0*, which is achieved exactly at the mechanical 
frequency u)2 ■ 

show ©/G against the frequency w in a logarithmic plot. 
As displacement sensitivities in current optomechanical 
devices are approaching 10“^®m/-\/Hz [TS1|T5], the appre¬ 
ciation of such small deviation of the quantum treatment 
from the classical effects appears to be not too far from 
being possible, albeit technically demanding. 

The small distance between SI and S2 will make it 
necessary to carefully study short-range and dispersion 
forces such as van der Waals interactions (vdW), where 
strength and power law scaling of the interaction with 
distance dramatically depend on the actual geometry of 
the two interacting systems [19] . Competing vdW effects 
are notorious for investigations of rather small gravity 
related effects and solutions for handling vdW have been 
already demonstrated [20] • 

Conclusions.- We have illustrated a method to infer the 
nature of the gravitational interaction between two mas¬ 
sive objects, in principle capable of discerning between 
a quantum and classical approach to gravity. Our ap¬ 
proach is based on the fundamental differences occurring 
in light of the possibility to prepare quantum coherent 
states of a system, within the quantum mechanical frame¬ 
work, which in turn gets manifested in the possibility to 
achieve coherent superpositions of distinguishable gravi¬ 
tational fields. Such a crucial difference between classical 
and quantum gravity can be revealed, in principle, in an 
optomechanical experiment, which showcases all the nec¬ 
essary ingredients to falsify one of the two treatments of 
gravity. 
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